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Abstract. Gravitational perturbations of the de Sitter spacetime are 
investigated using the Regge-Wheeler formalism. The set of perturbation 
equations is reduced to a single second order differential equation of the Heun- 
type for both electric and magnetic multipoles. The solution so obtained is used 
t' | , to study the deviation from an initially radial geodesic due to the perturbation. 

j^JTJ' The spectral properties of the perturbed metric are also analyzed. Finally, gauge- 

and tetrad-invariant first-order massless perturbations of any spin are explored 
following the approach of Tcukolsky. The existence of closed-form, i.e. Liouvillian, 
solutions to the radial part of the Teukolsky master equation is discussed. 

> 

o 

cn 

£f~) ' 1. Introduction 

o 

de Sitter spacetime is extensively investigated in the literature from different points 
of view. In fact, it is particularly relevant for inflation [T] and hence the physics of 
the very early universe. Furthermore, there are many recent experimental evidences 
and astronomical observations which indicate that the cosmological constant in the 
universe is positive leading to an accelerating expansion (see, e.g., Ref. [2] and 
references therein). The pioneering work of Gibbons and Hawking [3] extending to 
cosmological horizons the concepts of thermodinamics and particle creation previously 
introduced for black hole event horizons opened the way to the study of the properties 
of quantum field theory on de Sitter space, especially the temperature and entropy 
(see, e.g., Refs. [4j [5j [7]), and on asymptotically de Sitter spaces containing black 
holes. In this respect, it proved useful to express the de Sitter metric in static form 
(see, e.g., Ref. 0). 

We are interested here in studying metric perturbations in the static region of the 
de Sitter spacetime between the origin and the cosmological horizon. We will then use 
the Regge-Wheeler [9] formalism to decompose the perturbations in tensor harmonics 
and Fourier transform it with respect to time, as customary. We will account for 
perturbations of both parity, and investigate the features of geodesic motion in the 
perturbed gravitational field. In particular, we will consider an initially radial geodesic 
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in the unperturbed spacetime and study how the motion deviates from maintaining 
purely radial as a result of the perturbation. Furthermore, we will also consider 
gauge- and tetrad-invariant first-order massless perturbations of any spin following the 
approach of Teukolsky [10) . and discuss the existence of closed- form, i.e. Liouvillian, 
solutions to the radial equation. Note that our results are consistent with those of Refs. 
[TTj and [12] concerning the perturbations of the Schwarzschild-de Sitter spacetime, 
which however focused on the stability of the de Sitter spacetime in the presence of 
the black hole and on the scattering properties of both the event and cosmological 
horizons. 

Section 2 summarizes what is known about radial geodesies in de Sitter spacetime, 
while metric perturbations for both electric and magnetic multipoles are discussed in 
section 3. The equations for perturbations to radial geodesies are integrated in section 
4, while section 5 contains the derivation of Liouvillian solutions to the radial part 
of the Teukolsky master equation for massless perturbations of any spin. Concluding 
remarks are presented in section 6, while relevant details are given in the appendices. 



2. Radial geodesies in de Sitter spacetime 

Consider de Sitter spacetime, whose line element written in spherical-like coordinates 
is given by 

ds 2 = N 2 dt 2 - N~ 2 dr 2 - r 2 (dd 2 + sin 2 6»d0 2 ), (2.1) 

where N = (1 — i? 2 r 2 ) 1 / 2 denotes the "lapse" function. The spacetime region which 
can be accessed is the ball of radius 1/H and center at the origin. Timelike radial 
geodesies are characterized by their 4- velocity 

U = T{d t + r ) (2.2) 

with 

so that r denotes the proper time and 
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C = ± — VE^P; (2.4) 

E is the Killing constant associated with the stationarity of the metric. It is convenient 
to introduce an orthonormal frame adapted to the static observer family (at rest with 
respect to the coordinates) with basis vector fields 

e n = N- 1 d t , e f = Nd r: e § = l/rdg, = l/(rsin0)fy, (2.5) 
with dual basis of 1-forms 

up = Ndt, uj f = ~N- 1 dr, u) S = -rde, ^ = -rsm8d<j). (2.6) 
With respect to this frame the geodesic 4-velocity vector field can be written as 

U = 7 (e 6 + ve f ), (2.7) 

where 

7 = IW, v = N~ 2 (. (2.8) 
The parametric equations of the geodesies corresponding to initial values 

t(0) = , r(0) = 0, 6»(0) = 0, 0(0) = (2.9) 
can be written explicitly as follows. 
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(i) E 2 > 1. 

In this case it is convenient to use the notation E — cosh a. We have then 
'tai)h( ^IJt ) + sin] i n ^ i ^ tan] \ ( -7 // r ) sinlm 



arctanh ( — ; I + arctanh ' 



costw 



cosh( 



sinha . 
r(r) = — — — smh(iir) 
H 



(2.10) 



(ii) £ 2 = 1. 

In this case the initial conditions (|2 .9[) imply the trivial solution of a particle at 
rest at the origin 

r = 0, t = r. (2.11) 

(hi) E 2 < 1. 

In this case it is convenient to use the notation E = cos a. However r(0) 
cannot be made to vanish but it should be always greater than a certain value 
< r* — sin a/£f <1/H. We have then 



t(r) = — { arctanh 
H 



1 + sin a 



■ tanh 



cos a 







"1 


(h: 


+ arctanh 











tanh I -Ht 

cos a V 2 



r(r) 



if 



■ cosh(i7r) = r* cosh(fZr) 



(2.12) 



With a slight modification one can obtain the solution corresponding to more general 
initial data x a (0) = Xq. For instance, for r(0) 7^ the solution (|2.1ip for E 2 = 1 
becomes 



r = roe 



1 JV 



For a later use we also introduce the Newman-Penrose 13 frame 



1 iV 2 / 1 

l=J^d t +d r , n= — ^ — d t -d r 

The non- vanishing spin coefficients are 

1 N 2 V2 

p = --, H = - — , a = -{3 = -—cot&, 7 
r lr 4r 

The Weyl scalars are all identically zero. 



V2r 



H 2 r 



sin 8 



(2.13) 



(2.14) 



(2.15) 



3. Metric perturbations 

Let us consider the perturbations 

9nv = g^v + h^u (3.1) 

of de Sitter spacetime. Since the background metric (|2.1[) is static and spherically 
symmetric, we decompose the metric perturbation h^ v in tensor harmonics and Fourier 
transform it with respect to time, as customary. We will use the Regge- Wheeler [9] 
gauge to simplify the form of the perturbation. 
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3.1. Electric multipoles 

The metric perturbation h^ vi for the electric multipoles, is given by 



I 1 1 — 



" N 2 H 


Hi 












1 








sym 










e 


sym 


sym 


r 2 K 







sym 


sym 


sym 


r 2 sin 2 6K _ 





-iu)t \r 
HQ, 



(3.2) 



where the symbol "sym" inherited from [9] indicates that the missing components 
of h^v should be found from the symmetry h^ v — h Ufi and the functions Yiq are 
normalized spherical harmonics (see e.g. |14| ) with azimuthal index m 
by 



0, defined 



Y lQ = - 



1 21 + 1 



P/(cos( 



(3.3) 



L 



For I > 2, the system of radial equations we have to solve is the following (here 
1(1 + 1)): 



= H 1 ' + —(W + K)- 



W 

= K' - — 
r 



IL 



2ujr 2 
since 

H Q = H 2 = W, 
together with the algebraic relation 
(L-2 



H 



2rH 2 

"N 2 
K 

Tn 2 



-Hi, 







-W - 



-{H 2 L + 2lu 2 )H 1 



rN 2 



1 + t 2 uj 2 



N 



,L 



K. 



(3.4) 



(3.5) 



(3.6) 



Let us introduce the dimensionlcss variables x = Hr, ft = to/H, and set Hi = \H\ 
Solving for W in the constraint equation (|3.6[) and substituting into the system (|3.4 
we get (from now on a prime denotes differentiation with respect to x) 



H X ' 



K' = - 



(2il 2 + 3L - A)x 
(l-x 2 )(L-2) 
[2x 2 {L + 2n 2 ] 



Hi 



2n[{n 2 



L - l)x 2 



L] 



{l-x 2 ) 2 (L-2) 



L ( L ~3. Hl 



x{L + 2n 2 ) 



K, 



-K, 



(3.7) 
(3.8) 



2(L-2)nx 2 " x (l-a; 2 )(L-2)" 

which can be solved in terms of Heun's functions (see Appendix C). 

Figure Q] shows the behaviour of the metric perturbation functions Hi , K and 
W in terms of the rescaled radial variable x for I = 2 and different values of fi. The 
perturbative regime is generally lost very quickly. High frequency modes instead 
oscillate with small amplitudes over a wide range of the radial coordinate. The 
oscillations become very high when approaching the de Sitter horizon in the case 
of H\ and W . The amplitude remains instead practically constant in the case of K. 
In fact, in the high frequency limit ft ^> 1 one finds 

2x , tt 2 , . 

0, (3.9) 



K 



1 



(l-a; 2 ) : 



de Sitter spacetime: effects of metric perturbations on geodesies 



5 




(a) (b) 

Figure 1. The behaviour of the metric perturbation functions Hi, K and W in 
the case of electric multipoles is shown in terms of the rescaled radial variable x 
for I = 2 and different values of il. Fig. (a) corresponds to the choice Q = 1 
and the initial conditions are Hi(xq) = 10 — 4 = K{xq), with xq = 10 — 2 , so 
that W(xq) ft: 10 -4 . Fig. (b) corresponds to Q = 40 and initial conditions 
H 1 (x ) = 10 -4 = K(x ), with x = 10 -2 , so that W(x ) fa 7 x 10 -5 . Larger 
values of I correspond to a further growth of the perturbation functions. 



or equivalently 

^£+tfK = 0, (3.10) 

having defined z = arctanhx. The solution is then 

K = c\ sin(J7 arctanh x) + C2 cos(0 arctanha;), (3-11) 

where c\ and ci are arbitrary constants. 

A real solution for the metric can be obtained by considering W and K and H\ 
as real. In this case the non- vanishing metric components are 

g tt =N 2 (l + WY l0 cos ut), 
9tr = HxYio sin ut, 
g rr = - N~ 2 (l - WYioCosut), 
gee = - r 2 (l - KY l0 cos ujt), 

94,4> = ~ r 2 sin 2 0(1 - KY l0 coscut), (3.12) 

so that 

ds 2 = N 2 {1 + WYioCOSut) (dt - HxN^Yiosmujtdry 

- N- 2 {1 - WY l0 coscut)dr 2 

- r 2 (l - KY m cos ut){d6 2 + sin 2 9d(j) 2 ), (3.13) 
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to first order in the perturbation quantities. A natural orthonormal frame associated 
with this form of the metric is then 

O 6 = N (l + ^WY m cos utj (dt + HiN~ 2 Yi sinutdr 

= N ( X + \ WYl ° COSCJ t) dt + H 1 N~ 1 Y m sinujtdr, 
Q? = - N- 1 (l - lwY lo cosuit] dr, 
Q § = 
= 



-r I 1 - -KY m cosujt | de. 
— rs'mO ( 1 — -KYiq cos _,_■/ j <!<>. 



(3.14) 



One can then introduce a Newman-Penrose |13) frame in a standard way [15] 

I = (tt 5 + fT f )/\/2, n= (fi -^)/V2, m= (O e ~ + i^)/V2, (3.15) 

and conveniently rotate it according to I — > ly/2/N, n — > nN/y/2. We list here the 
non- vanishing first-order spin coefficients and Weyl scalars (hereafter C = coswi and 
S = sintut): 

r lr 
N 2 N 2 



P 



Z + ~{rK' - W)Y l0 C - ^KY l0 S, 



-(rK'-W)Y l0 C+-KY l0 S, 
2r 4r 4 



V2 



V2„ fdYi 



, cot6> + ^— K 
4r 8r 

V2~dY l0 
p = — a — Hi — a, 



de 



-ro\HY„, )C-^HxY l0 S, 
8r 



7 = 



4r d6 

_ l (- N * W ' + 2uR x + 2rH 2 W)Y lQ C - ^WY l0 S, 



v = ^ { WC-H lS ) dY " 



8r 
Yo 



de ' 



AN 2 



l(N 2 W' - 2uH l )C + ljWS], 



and 



■04 



2N 2 r 
1 



2N 2 i 



de ' 

(WC + H 1 S) 2 Y°i, 



(3.16) 



J2_ 

AN 2 r 2 

L 



(HiLN 2 + 2ruK) 
2~ruj 



C + (H X N 2 + rujK)S 



Aujr c 



{HtN' + rwIQYmC, 



V2 
8r 2 

8T 2 



(HiLN 2 + 2ruoK) 
2ru) 



(wc-ms^Y^, 



C - (HiN 2 + ru;K)S 



dYio 

de ' 



dYio 

de ' 



(3.17) 
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X 



Figure 2. The behaviour of the spectral index S in the case of electric multipoles 
is shown in terms of the rescaled radial variable x for fixed values of t = 7r/4, 
8 = n/2, 1 = 2 and Q = 1. The choice of initial conditions is the same as in 
Fig. [T] (a). The variation of radial coordinate is limited to the range where the 
perturbative regime holds. 



where 



2^ cote + LY l0 . 
du 



(3.18) 



Note that the first-order Ricci scalar coincides with the zeroth-order one. 

It is useful to introduce the speciality index S, which is a suitable combination of 
Weyl scalars defined by 



5 _ 27J 2 _ 27 (^o^4-^ 2 3 ) 2 



(^4 + 3^ 2 2 ) 3 ' 



(3.19) 



and has the value 1 for algebraically special spacetimes. Note that the imaginary part 
of S turns out to be identically vanishing, so that it is a purely real quantity. Its 
behaviour as a function of the rescaled radial coordinate x for fixed values of t and 6 
is shown in Fig. [5J 



3.2. Magnetic multipoles 

The metric perturbation h^, for magnetic multipoles, is given by 












ho 


sym 








In 


sym 


sym 








sym 


sym 


sym 






e~ iut sin 



} dY l0 

de ' 



For I > 2, the system of radial equations we have to solve is 



= h ' - 



N 2 (L - 2) 



, 2rH 2 iuj 



(3.20) 



(3.21) 
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(a) 



(b) 



Figure 3. The behaviour of the metric perturbation functions ho and hi in 
the case of magnetic multipoles is shown in terms of the rescaled radial variable 
x for I = 2 and different values of Q. The initial conditions are chosen so that 



f! = 1 and Q 



: = hi(xo), with xo 
40, respectively. 



iir 



Figures (a) and (b) correspond to 



where again L = 1(1 + 1). 

By introducing the dimensionless variables x = Hr, f2 = lu/H as above and 
setting ho = iho, the system p.21| becomes 



- , 2 - [{n 2 + L- 2)x 2 +2-L] 
ha — —n-a 2 "ii 



v = 



7 h 



'Ii- 



l—a: 



-M, 



(3.22) 
(3.23) 



(l-z 2 ) 2 

which can be solved in terms of Heun's functions (see appendix C). 

Figure [3] shows the behaviour of the metric perturbation functions ho and h\ in 
terms of the rescaled radial variable x for I = 2 and different values of Vt. The same 
considerations holding for electric multipoles apply also in this case. 

A real solution for the metric can be obtained by considering ho and h\ as real. 
In this case the non- vanishing metric components are 



at2 l ■ . ■ /i^O , , . n dYio 

gtt = N , gt<p = h smujt smO— — , g r<j> = hi cosut sm8— — . 

do do 



9rr = — N 

so that 



94>4> 



-r 2 sin 2 i 



ds 2 = ds 2 



li i S] -t- 2 sin 9 ^° (jiQ sin utdcf) + 2h\ cos wttir^ 



(3.24) 



(3.25) 



to first order in the perturbation quantities. A natural orthonormal frame associated 
with this form of the metric is then 

fl° = Ndt + N^ 1 smO^-h smLjtd<j), 
dO 



de Sitter spacetime: effects of metric perturbations on geodesies 



ff = -N~ 1 dr + Nsiae^-hicosLJtdd), 

dO 

n § = - rd9, = -rsi 



(3.26) 

One can then introduce a Newman-Penrose frame in a standard way, as in Eq. (|3. 1 5[) . 
We find that the non-vanishing first-order spin coefficients, Weyl scalars and Ricci 
terms are as follows: 



f 



L 



[iV 2 /iiC + h S]Y l0 , 



P r ' 2 N 2 r 2 

N 2 T 

M= - — +i T ^[N 2 h i C-h S]Y l0 , 
2r 4r z 



V2 .V2 

——COtO - l—r 

4r 4r 2 



2ujrN 2 



dY m 

de 



. V2N 2 dY l0 



7 



^ +^[^^0 -h S]Y l0 , 



/2 

i— — I \N 2 H 2 rh\ + ujho\C 
f 



2ur 2 
L 



AN 2 r 2 



[2ujrh + N 2 {N 2 (L - 2) - 2u 2 r 2 )h 1 ]S 



dY l0 

de 



. V2_ 



{[N 



2 H 2 rh 1 +uh ]C 



+ ? [Zurho + N 2 {N 2 (L - 2) - 2u 2 r 2 )h 1 ]S 



2bJT 



dYm 

de ' 



and 



1p2 



H 2 rhi +u>h ]C 



+ 



2ur 2 



[2wrh + N 2 (N 2 (L - 2) - 2uj 2 r 2 )h 1 ]S \ 2 Y°i 



.V2(L - 2) 



4r 3 
N 2 
~4ujr 4 



hiC- 



rN 2 



[N 4 hi-Lorh }S 



dY, 



/o 



de 



i L{L-2)h l Y l0 S, 



i>3 = i 



1p4 



.V2N 2 (L-2) 



8r 3 



rN 2 



[7V 4 /n - urhojS 



dY l0 

de ' 



H 2 rhi +uh ]C 



■—r[2wrh Q + N 2 (N 2 (L - 2) - 2w 2 r 2 )h 1 }S \ 2 Y°i. 



(3.27) 



(3.28) 



The behaviour of the speciality index (|3.f 9[) as a function of the rescaled radial 
coordinate x for fixed values of t and e is shown in Fig. 01 Note that S is purely real 
also in this case. 
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Figure 4. The behaviour of the spectral index S in the case of magnetic 
multipoles is shown in terms of the rescaled radial variable x for fixed values 
of t = 7r/4, 8 = 7r/2, I = 2 and Q = 1. The choice of initial conditions is the same 
as in Fig. [3] (a). The variation of radial coordinate is limited to the range where 
the perturbative regime holds. 



4. Perturbations of radial geodesies 

We are now in a position to consider the perturbations 

U = U gco + u (4.1) 

to the radial geodesic (|2.2|) . We are interested in finding the solution for u such that 
the path is geodesic in the perturbed gravitational field, i.e. VjyU = 0, to first order 
in the perturbation. 



4-1- Electric multipoles 

The perturbed 4-velocity is given by 



-^Y lo d t +CxN 2 Y lo d r + %^d ( 



N 



d9 



where T\, d and £i are functions of r. Real solutions for u require 



£i = A\ cos Lot + Ai sinwt, 



Mo 



N 2 , 

B\ H — —W ) cosw£ + B 2 smut 



l r 

E 



(4.2) 



(B 1 + E 2 W) cosout + [B 2 + ETqCqHA smut (4.3) 



where the functions A\ , A2 and B\ , B2 satisfy the systems of equations 



dAi 
dr 



2E 2 - N 2 , 



2Euj 



-W 



dA 2 uj . , 



(4.4) 
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and 



dr 
dB 2 
dr 



Co" 



^B 2 + £ToCo-ffi 



2E 2 - N 2 



-W 



(4.5) 



respectively, which can be integrated numerically. 

The behavior of the perturbed velocity functions (|4.3[) is shown in Fig. [5] in terms 
of the rescaled radial variable x for different values of t and fixed values of I, and 
E. 



4-2. Magnetic multipoles 

The perturbed 4- velocity is given by 

6 dY l0 



u 



-iu)t 



r 2 sin 9 dO 

where £i is a function of r. Real solutions for u require 
where the functions , A 2 satisfy the system of equations 



dAi 
dr 

dA 2 
dr 



Co 



A 2 



2E 2 - N 2 

'en 2 



ho 



rH 2 N' 



hi 



to 
Co 



2r Co 

uirN 2 



h - uirN 2 1 



N 2 (L 



2u 2 r 



hi 



(4.6) 
(4.7) 



(4.8) 



which can be integrated numerically. 

The behavior of the perturbed velocity function (14. 7[) is shown in Fig. [5] in terms 
of the rescaled radial variable x for different values of t and fixed values of I, fi and 
E. 



5. Teukolsky equation 



In this section we study gauge- and tetrad-invariant first-order massless perturbations 
of any spin (s = ±2 for gravitational waves, s = ±1 for electromagnetic waves, 
and s = ±1/2 for neutrinos) in the de Sitter background following the approach of 
Teukolsky [10]. Although this problem have been already addressed in the literature 
(see, e.g., Ref. [16], where the perturbation equations for fields of different spin 
have been solved by using the technique of Dcbye potentials), we are are interested 
here in investigating the existence of closed-form, i.e. Liouvillian, solutions to the 
perturbation equations. 

The Teukolsky master equation, in this case, reads as 

N 2 







/J- 

\27V 2 



d t , 



-d r 



1 



1 



2r 2 
-2N 2 ) 



2r 2 sin 2 9 
cot 



0, b 



rN 2 
is cos 6 
r 2 sin 2 i 



1 

2? 1 



' 2r 2 
(l + s)[3s + 2-27V 2 (2s 



1)] 



sm 



(5.1) 
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(c) 



Figure 5. The behavior of the perturbed velocity functions e -1 "'^!, e — 
and e — IW *Fi in the case of electric multipoles is shown in terms of the rescaled 
radial variable x for different values of t and the same choice of parameters as 
well as initial conditions for the perturbed metric functions as in Fig. [T](a). The 
functions A\,A^ and B\,B2 are taken instead identically vanishing at the initial 
point. The unperturbed geodesic radial motion is directed outwards with energy 
parameter E = 1.1. Note that in the high frequency limit the perturbed velocity 
functions exhibit instead an oscillating behavior similar to that of the perturbed 
metric functions, the amplitude of c" 1 "'^ at fixed distance being about 10 2 
smaller than that associated with f i and Ti . 
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Figure 6. The behavior of the perturbed velocity function e -Ki,t £i in the case 
of magnetic multipoles is shown in terms of the rescaled radial variable x for 
different values of t and the same choice of parameters as well as initial conditions 
for the perturbed metric functions as in Fig. [3] (a). The functions A\,A2 are 
taken instead identically vanishing at the initial point. The unperturbed geodesic 
radial motion is directed outwards with energy parameter E = 1.1. Note that 
in the high frequency limit the perturbed velocity function exhibits instead an 
oscillating behavior similar to that of the perturbed metric functions. 



Note that this equation is solved by ip = ipo with s = 2 and by ip = ipi/ P with ,s = —2, 
where ipo and -04 are the perturbed Weyl scalars given by Eqs. (|3.17[) and (|3.28j) for 
electric and magnetic multipoles, respectively. Therefore, the latter equations provide 
a link between the perturbed metric functions and the solutions to the Teukolsky 
equation. 

Let us now study the Teukolsky equation in completely general form, relaxing 
the assumption of azimuthal symmetry. The master equation (15 . 1[) admits separable 
solutions of the form 



-\u)t imcp 



R(r)S(8), 



(5.2) 



where uj > is the wave frequency and m is the azimuthal separation constant. 

The angular function satisfies the equation of the spin-weighted spherical 
harmonics, i.e. 

Id/. ,dS(0)\ 



sin e de 



sin t 



d9 



V, 



where 



(ang) 



(9) = L - s 2 - 



(s COS 



i;mg) (0)S(0)=O, 
m) 2 



(5.3) 



sm 



and the separation constant has been chosen according to Teukolsky. 
The radial equation can be cast in the form 



dr \ ar 



+ V {lad) (r)R(r) = 0, 



(5.4) 



(5.5) 



where Q 



r 2 N 2 and 

1 f , 2iws N 2 
N 4 



N 4 



[L + (1 + s)(3s + 2)] + 2 — (1 + fl )(2a + 1) } . (5.6) 
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5.1. Liouvillian solutions to the radial equation 



The general solution to the radial part of the Teukolsky equation can be expressed 
in terms of Heun functions. We are interested here in investigating the existence 
of closed-form, i.e. Liouvillian, solutions to the radial equation (15 ,5|) governing the 
perturbations. Roughly speaking, the set of Liouvillian functions includes the usual 
elementary functions such as exponential, trigonometric, logarithmic functions, etc., 
but not generic hypergeometric functions or other special functions. We apply the 
well-known Kovacic algorithm |17) outlined in appendix D. 

The radial equation (15.51) can be converted in normal form, i.e. 

<1 "' /( " -W {Tad) (r)y(r), 



dr 2 

by the scaling 

R = yQ- {1+s)/ \ 
where the function W( ra( j) is given by 



W(wd)(r) = 
with parameters 
c = L, 
do = — 2iuis 



H 



r 2 (r - r H ) 2 



H 



(r + r H ) 2 



r 



'H 



'H 



d% 



T-. — [w 2 r^ =F 4iwrH-s 



(r-r H ) {r + r H y 



s 2 - 2L+1}. 



(5.7) 
(5.8) 

(5.9) 



(5.10) 



Let us discuss in detail the case of type 1 solutions (see again appendix D). For all 
possible pole structures of (15.91) . all functions fl generated by the Kovacic algorithm 
have the form 



Oo 

r 



<H 



(r + r H ) (r-rii)' 



where 



«o=2 [1 + <5o\/l +4c J , 



1 



'H 



l + 4c 



H 



(5.11) 



(5.12) 



the quantities So, 5^ and djj taking independently values ±1. 

The necessary conditions for type 1 solutions are fulfilled if the following 
constraint is imposed: 



d 



0. 



implying that the order of the pole at infinity is 2. 
The algorithmic condition d — n is given by 



ao — ctjj — a H = n, 



where 



1 + See^l + 4(c + c+ +c H ) 



±1. 



(5.13) 



(5.14) 



(5.15) 



For every pole structure, we construct the Liouvillian solutions to Eq. (I5.7[) by 
deriving a consistent recursion relation which generates the finite set of coefficients in 
the polynomial V . To obtain a recursion relation with the fewest number of terms 
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we take, without loss of generality, V to be expanded about a pole of non-vanishing 
order, say r = 0, as follows: 



a k r 



a n ^ 0. 



(5.16) 



k=0 



After substituting this expression for V into Eq. (|D.2[) we obtain in the general case, 
when all three poles are present, the recursion relation 







(k - l)(k - 2 - 2n + 2 aoo ) - r H (g+ - g H ) 



2k(a+ - a H ) - r H (g+ + g H ) 



a k 



where 



90 



9% 



- (k + l)r 2 H (k + 2a Q )a k +i, 



d + 2 — (aj { - a H ), 
d%T—(2a + a T), 

iff 



fc = 0, . 



with the understanding that a; = if i < —1 or i > n + 1. 
Liouvillian solutions (ID.3|) to Eq. (|5.7[) are thus given by 



y = Vr ao (r + r H ) a H (r H - r) a H , 



implying that 



R = Vr 1 + s r aa 



2 (r + rjy) H 2 (rH~r) H 2 



(5.17) 



(5.18) 



(5.19) 



(5.20) 



6. Concluding remarks 

The present work has been motivated by a previous study on gravitational waves 
about de Sitter backgrounds [IB], which is relevant for the stochastic background 
emerging from inflation, and by the analysis in [19] of gravitational waves and their 
geodesies. We have solved the gravitational perturbation problem in the static region 
of the de Sitter spacetime between the origin and the cosmological horizon by using 
the Regge-Wheeler formalism. The set of perturbation equations has been reduced 
to a single second order differential equation of the Heun-type for both electric and 
magnetic multipoles. We have then investigated the geodesic motion in the perturbed 
gravitational field. In particular, an initially radial geodesic in the unperturbed 
spacetime cannot remain radial as a result of the perturbation. The deviation from 
radial motion has been shown to occur in both polar and azimuthal directions. Finally, 
we have investigated gauge- and tetrad-invariant first-order massless perturbations of 
any spin following the approach of Teukolsky. The general solution to the radial part 
of the Teukolsky equation can be expressed in terms of Heun functions. However, we 
have shown that there also exists in addition a class of closed-form, i.e. Liouvillian, 
solutions. 
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Appendix A. Electric multipoles: / = 0, 1 
Appendix A.l. I = 

The system of perturbation equations in the case I = is given by 



2w ~ 



rH< 



= K' 



#2 
r 



K 

rN 2 ' 



We can set 

H Q = H 2 = W, Hi = 0, 
without loss of generality, so that Eq. (|A.1|) reduces to 

rH 2 



= W' 



= K' 



N 2-W ~ + H 2 )K, 



W 



K 



r rN 2 

In terms of the dimensionless variable x = Hr the corresponding solution is 



K = 



W = 



ji^x 2 



^fT+W 



Cl 



1 



x 



l + x 



- C '2 



VT^X 2 

where fl = ui/H. 



-ci 



1 - x 



Vl+0 2 /2 



+ c 2 



1-x 
l + x 

1 - 



-Vl+fi 2 /2' 



1 



/2' 



Appendix A. 2. 1 = 1 

The system of perturbation equations in the case Z = 1 is given by 
= JTo' - 



rTV 2 
H2 

1 

r 



N 4 N 2 ' 



= K' 

plus the algebraic relation 
tor 

* - w K - 

We can set 

H = H 2 = W, 
so that Eq. (|A.5j) becomes 

= a" - ^. 

r 

In terms of the dimensionless variable x = Hr the corresponding solution is 
K = c\ sin(i7 arctanha;) + c 2 cos(f2 arctanha;), 



(1 + H 2 r 2 
rN 2 



■W, 



w = 



fix 



(1-x 2 ) 



[ci cos(f2 arctanha;) — c 2 sin(51 arctanha;)] 
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Appendix B. Magnetic multipoles: / = 1 

The system of perturbation equations in the case I = 1 reduces to 

= h '- — +uhi. (B.l) 
r 

Without loss of generality we can set hi = 0, so that the previous equation becomes 

= ^'-—, (B.2) 
r 

whose solution is simply 

h = cir 2 . (B.3) 

Appendix C. Decoupling the perturbation equations 

The system of perturbation equations of both parity can be reduced to a second- 
order differential equation for a single perturbation function, the remaining ones being 
determined either algebraically or by differentiation. Such a master equation can be 
cast in general in the form of a Heun equation [501 HI] 

d 2 y | / 7 | S | 6 \ dy afi-q y — Q (CI) 

dz 2 \z z — 1 z — a j dz z(z — l)(z — a) 

where {a, /3, 7, 5, e, a, q} are complex arbitrary parameters such that j+S+e = a+/3+l 
and a ^ 0,1. This equation has four regular singular points, at {0, 1, a, 00}. 

This is a generalization of the hypergeometric equation. In fact, in the special 
case where e = and q — a/3a, Eq. (jC.ip becomes 

>S + [(a + /> + l)-7$ 
dz z dz 

which is just the hypergeometric equation in its standard form, with solution 

V = ciF(a, 0; 7; z) + c 2 z 1 -^F{a - 7 + 1, £ - 7 + 1; 2 - 7; z). (C.3) 

Appendix C.l. Electric multipoles 

The coupled system p.7p ~ (|3.8p is equivalent to solve a second-order differential 
equation for K of the Heun type. The procedure is as follows. First, solve algebraically 
Eq. (13.8[) for Hi and substitute it into Eq. (|3.7p . thus obtaining a second-order 
differential equation for K . Introduce then the new variable z = x 2 togheter with the 
rescaling 

K(z) = f(z)y(z), f(z)=z a Hz-ir(z-ar, a =~ -^^ . (C.4) 

The function y thus satisfies an equation of the form (jC.lj) with parameters 

3 

7 = 2o?i +2' ^ = ^ a2 1' e ~ — 1' 
a = ai+a 2 + a 3 + --—, p = a + —, 

q = 3(ori + a 2 + 1)(1 - a) + a(a/3 - 2) - 3a 2 + -, (C.5) 



*(* " !)tt + [(« + Z 3 + 1) - 7)37 + a Pv = 0. ( C - 2 ) 
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and 

a 1 = ~ + ^VT+4L, a 2 = -ia 2 ^, a 3 = 2, (C.6) 

with vT+4L = 21 + 1 and a±, a 2 , 03 taking independently the values ±1. The signs 
have to be chosen so as to avoid singular behaviours at z = and z = 1. 



Appendix C.2. Magnetic multipoles 

The same procedure outlined above for the case of electric multipoles can be followed 
here to reduce the coupled system (|3.22l) ~ (|3.23[) to a second-order differential equation 
for h\ of the Heun type. First, eliminate ho in favour of hi, then introduce the new 
variable z = x 2 and a rescaling h\(z) = f(z)y(z) of the form (|C.4[) . The function y 
thus satisfies an equation of the form (|C.1|) with parameters 

7 = 2ai — -, 5 = 2a 2 +3, e = 2a 3 , 
a = ai+a 2 + --— : f) = a+—, 

q = apa, (C.7) 

and 

3 cr £1 

ai = - + -^vTTiZ, a 2 = -l-io-2-, a 3 = 0, (C.8) 

implying that e = 0. Therefore in this case the Heun equation reduces to the 
hypergeometric equation (|C.2|) . 



Appendix D. The Kovacic algorithm: an overview 

The Kovacic algorithm applies to a second-order ordinary differential equation of 
normal form 

S = W{r)y, (D.l) 
dr z 

when W(r) is a given element of C(r), the field of rational functions with coefficients 
in the field of complex numbers. Kovacic proved that all Liouvillian solutions to 
Eq. (|D.lj) belong to three mutually exclusive types. Furthermore, for each type, he 
provided an algorithm which determines whether a Liouvillian solution exists for that 
type, and constructs the solution in case it does exist. The types are as follows: 

Type 1. The differential equation has a solution of the form e^ sldr , where tl € C(r). 

Type 2. The differential equation has a solution of the form e^ odr , where Cl is algebraic 
over C(r) of degree 2 (i.e., Cl is a solution of a polynomial equation of degree 2 
with coefficients in C(r)), and solutions of type 1 form an empty set. 

Type 3. All solutions of the differential equation are algebraic over C(r) of degree 2, and 
solutions of type 1 and 2 form two empty sets. 

If none of the previous cases applies, the differential equation has no Liouvillian 
solution. The construction of each type of Liouvillian solution is derived from the 
pole structure of the rational function. Kovacic also discussed some conditions that 
are necessary for cases 1, 2 or 3 to hold. For type 1 every pole of W(r) should have 
even order or else order 1, and the pole at 00 should have either even order or order 
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greater than 2. For type 2 a necessary condition is that W(r) must have at least one 
pole that has either odd order greater than 2 or order 2. Last, for type 3 the order of 
a pole of W(r) cannot exceed 2, and the order of the pole at oo must be at least 2. 

Let us discuss in detail the case 1. The remaining cases 2 and 3 apply as well and 
can be treated similarly. 

For a fixed rational function W(r), type 1 Liouvillian solutions of Eq. (|D.lj) are 
determined in the following way. A finite set of constants d is derived algorithmically 
from the pole structure of W(r) and for each d a function f2(r) is a l so constructed. 
If any d is a non-negative integer n and, for the corresponding f2(r), the differential 
equation 



dr 2 dr 



^ + n 2 - W(r) 
dr 



T{r) = (D.2) 



has a polynomial solution V of degree n, then Eq. (|D.1|) has a Liouvillian solution 
given by 

y = Vef hdr . (D.3) 

If this is not the case for any d, then Eq. (jD.ll) has no type 1 solutions. 

Since W(r) depends on different independent parameters, several different pole 
structures may occur by imposing conditions on the parameters which alter the order 
of poles. For each different pole structure, in general, d and fl are expressed in terms 
of the parameters. Solutions to Eq. (|D.1[) may be obtained by requiring d = n, which 
then is just a constraint on the parameters, and finding all sets of parameters for which 
Eq. (|D.2|) has a polynomial solution of degree n. 
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